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Abstract 


We describe the isolation, biological and genetic characterization of a host-range variant of bovine coronavirus (BCoV) detected in water 
buffalo (Bubalus bubalis). By conventional and real-time RT-PCR assays, the virus was demonstrated in the intestinal contents of two 20-day-old 
buffalo calves dead of a severe form of enteritis and in the feces of additional 17 buffalo calves with diarrhea. Virus isolation, hemagglutination 
and receptor-destroying enzyme activity showed that the buffalo coronavirus (BuCoV) is closely related to BCoV but possesses some different 
biological properties. Sequence and phylogenetic analyses of the 3’ end (9.6 kb) of the BuCoV RNA revealed a genomic organization typical of 
group 2 coronaviruses. Moreover, the genetic distance between BuCoV and BCoV was proven to be the same or even higher than the distance 
between other ruminant coronaviruses and BCoV. In conclusion, our data support the existence of a host-range variant of BCoV associated with 


enteritis in buffaloes. 
© 2007 Elsevier Inc. All rights reserved. 


Keywords: Buffalo; Enteritis; Bubaline coronavirus; Biological characterization; Molecular characterization 


Introduction 


Coronaviruses (CoVs) (order Nidovirales, family Coronavir- 
idae) are enveloped, positive-sense single-stranded RNA particles 
that are responsible for enteric and/or respiratory disease in 
mammals and birds. Their ability to evolve through genetic 
recombination and/or point mutation is well-known so that they 
give rise to new viral genotypes or mutants with different tissue or 
host tropism. Currently, CoVs are organized into three antigenic 
groups with group 2 including bovine-like (subgroup 2a) and 
SARS-like (subgroup 2b) viruses. Bovine coronavirus (BCoV) 
belongs to subgroup 2a together with mouse hepatitis virus, 
sialodacryadenitis virus, porcine hemagglutinating encephalo- 
myelitis virus, human coronavirus (HCoV) OC43, human enteric 
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coronavirus (HECV) 4408 (Enjuanes et al., 2000) and the newly 
recognized equine coronavirus (ECoV) (Guy et al., 2000), HCoV- 
HKU1 (Woo et al., 2005), and canine respiratory coronavirus 
(CRCoV) (Decaro et al., 2007a; Erles et al., 2003). There are 
multiple genetic and antigenic evidences that several subgroup 2a 
CoVs, such as HCoV-OC43, HECV-4408, PHEV and CRCoV, 
have arisen as a consequence of trans-species infections caused by 
BCoV (Zhang et al., 1994; Vijgen et al., 2005, 2006; Erles et al., 
2007). BCoV can cause severe diarrhea in newborn calves 
(Snodgrass et al., 1986), “winter dysentery” in adult cows (Cho 
et al., 2000; Saif et al., 1991) and respiratory tract illness in calves 
and cows (Lathrop et al., 2000; Storz et al., 2000). The same virus 
strain could be responsible for simultaneous appearance of enteric 
and respiratory disease in the same animals (Chouljenko et al., 
2001) as well as in both calves and cows (Travén et al., 2001). 
BCoV infection is responsible for severe economic losses in cattle 
herds due to decreased milk production in dairy cows and reduced 
weight gain in beef calves (Saif et al., 1998), with generally low 
mortality rates. 
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Water buffalo (Bubalus bubalis) represents a ruminant spe- 
cies important in the economy of several countries, including 
Brazil, India, Vietnam and some regions of southern Italy, 
where the breed Mediterranean Italian buffalo produces high 
quality milk employed for production of the buffalo “mozza- 
rella cheese” (Romano et al., 2001). Recently, bovine-like CoVs 
have been identified in wild or domesticated ruminants, 
including several species of deer, waterbuck antelope (Tsune- 
mitsu et al., 1995), giraffe (Giraffa camelopardalis) (Hasoksuz 
et al., 2007), alpaca (Lama pacos) (Jin et al., 2007) and sable 
antelope (Hippotragus niger) (Spiro D., R. Halpin, S. Wang, M. 
Hasoksuz, X. Zhang, K. Alekseev, A. Vlasova, D. Janies, E. 
Ghedin, and L. Saif, unpublished). However, CoVs have never 
been isolated from buffaloes, although there is a single report 
on the detection of BCoV antibodies in Bulgarian buffaloes 
by hemagglutination inhibition and virus neutralization tests 
(Muniiappa et al., 1985). In the same report, hemagglutinating 
activity which could be suppressed by a BCoV-specific bovine 
serum was also demonstrated in some bubaline fecal samples. In 
this paper, we describe the isolation and biological and genomic 
characterization of a bovine-like coronavirus from a buffalo 
herd affected by severe diarrhea and calf mortality. 


Results 
Clinical outbreak 


The outbreak occurred between October 2006 and April 
2007 in a herd of Mediterranean Italian buffalo (B. bubalis) in 
Campania (southern Italy). No cattle, sheep or goat farms were 
present in the vicinity of the outbreak. At that time the herd 
consisted of 460 buffaloes, including 215 lactating cows, all 
vaccinated against colibacillosis, clostridiosis and salmonello- 
sis. Buffalo calves were removed from their dams shortly after 


Table 1 
Bovine-like CoV RNA titers detected by real-time RT-PCR in buffalo calves 
with diarrhea 


Prot. no. Age of the Sample RNA copies/ul of 
animals type template 
179/07-A 20 days Intestinal content 2.13108 
179/07-B 20 days Intestinal content 9.41 x 10° 
179/07-1 8 days Feces 1.25x 10° 
179/07-2 2 months Feces 4.10x10* 
179/07-3 2 months Feces 5.98 x 10° 
179/07-4 5 days Feces 4.55x 104 
179/07-5 1 month Feces 8.14« 104 
179/07-6 2 months Feces 3.27x 10° 
179/07-7 5 days Feces 2.90 x 107 
179/07-8 12 days Feces 6.71 10° 
179/07-9 2 months Feces 2.76 x 10° 
179/07-10 2 months Feces 7.18 x 107 
179/07-11 3 months Feces 5.23x 107 
179/07-12 1 month Feces 2.04 x 10° 
179/07-13 1 month Feces 1.77x 10° 
179/07-14 2 months Feces 3.34 10" 
179/07-15 1.5 months Feces 4.29 x 107 
179/07-16 2 months Feces 6.83 x 10° 
179/07-17 2 months Feces 2.69 x 104 


Fig. 1. Virus isolation of the bubaline CoV on human rectal tumor cells. 
(A) Cytopathic effect (syncytia) caused by CoV strain 179/07-11 (72 h 
postinoculation). (B) Cytoplasmic fluorescence detected by the immunofluores- 
cence assay using a BCoV-specific serum (24 h postinoculation). 


their birth and placed in separate hutches according to the 
gender where they were hand-fed fresh colustrum for 5 days. 
Neonatal mortality was firstly observed in October 2006 in 30 
5- to 20-day-old calves (out of 40 newborns) that displayed 
severe diarrhea and died despite treatment with antibiotics 
(oxytetracyclin and amoxicillin). Simultaneously, gastroenteric 
disease was also observed in older calves (1-3 months of age). 
On April 2007, neonatal mortality and enteric signs in calves 
persisted and two carcasses of 20-day-old dead calves together 
with fecal samples from additional 17 diseased calves were 
submitted to our laboratory for routine analysis. 


Identification of a bovine-like CoV in buffalo calves 


At necropsy, the carcasses of the two dead buffalo calves 
showed severe gastroenteritis, with enlargement of the mesen- 
teric lymph nodes and gall bladder. By conventional RT-PCR 
(Erles et al., 2003), bovine-like CoV RNA was detected in the 
intestinal content of the dead animals as well as in all 17 fecal 
samples from calves with diarrhea. Sequence analysis of the S- 
gene fragments showed a 100% nucleotide identity among the 
CoV strains detected in different animals. 
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Table 2 
Hemagglutination (HA) and receptor-destroying enzyme (RDE) titers of bubaline CoV 179/07-11 and BCoV strains 
CoV strain Origin HA titer® RDE titer* 

4°C SP °C 

Mouse Chicken Mouse Chicken Mouse Chicken 
BuCoV-179/07-11 Feces 512 - 128 - 128 - 
BCoV-339/06 Feces 128 64 64 32 64 16 
BCoV-438/06-2 Nasal swab 512 256 128 16 128 64 
BCoV-9WBL7 Unknown 1024 256 1024 64 256 = 


* HA and RDE titers per 50 11 of fecal suspension are expressed as reciprocals of the highest dilutions producing HA and resulting in complete disappearance of HA, 


respectively. Titers <2 are indicated as negative (—). 


Gastrointestinal parasites of ruminants, including Cryptos- 
poridium spp., were not detected by zinc sulfate flotation or 
Ziehl Nielsen staining, whereas by bacteriological investigations 
only common bacteria with poor pathogenic significance 
(E. coli, Streptococcus spp., Staphylococcus spp., Corynebac- 
terium spp.) were isolated from the intestinal contents or fecal 
samples. Using real-time RT-PCR, bovine-like CoV RNA was 
detected at low titers in the intestinal contents of the dead 
animals, whereas higher viral loads were found in the fecal 
samples of calves with diarrhea, with a peak of 5.23 x 10’ RNA 
copies/l of template in calf 179/07-11 (Table 1). 


Biological characterization of the bubaline CoV 


The fecal sample of calf 179/07-11 was used for virus isolation 
attempts on human rectal tumor (HRT-18) and Madin Darby 
bovine kidney (MDBK) cells. Cytopathic effect (CPE) consisting 
of syncytia formation and subsequent cell lysis was evident at the 
2nd passage on HRT-18 cells (Fig. 1A), where the growth of a 
bovine-like CoV was confirmed by the cytoplasmic fluorescence 
detected by the IF assay using a BCoV-specific serum (Fig. 1B). 
At the Sth passage on HRT cells the viral titer reached values of 
10°? TCIDs9/50 pl of viral suspension and 1.97 10° RNA 
copies/ml of template, as calculated by the Karber’s method and 
real-time RT-PCR, respectively. Replication of strain 179/07-11 
in MDBK cells was not as efficient as in HRT-18 cells, 
considering that CPE was not observed on MDBK monolayers 
and only few positive cells were detected by the IF assay at the | st 
passage, with viral titers (as calculated by real-time RT-PCR) 
decreasing progressively in the subsequent passages until defi- 


RdRp HE 


nitive loss of the replication ability. MDBK cells were confirmed 
to be poorly permissive to 179/07-11 replication even when the 
HRT-18 cell-cultured virus was inoculated. 

The isolated strain was evaluated by assessment of hemag- 
glutination (HA) and receptor-destroying enzyme (RDE) 
activities in comparison with the BCoV isolates 339/06 (Decaro 
et al., 2007b), 438/06-2 (Decaro et al., in press) and 9)WBL7 
(kindly provided by Dr Paolo Cordioli, Istituto Zooprofilattico 
Sperimentale della Lombardia e dell’ Emilia Romagna, Brescia, 
Italy). All isolates agglutinated mouse erythrocytes at titers 
higher at 4 °C than at 37 °C. The bovine isolates caused HA with 
chicken erythrocytes at low titers at both 4 °C and 37 °C, 
whereas the bubaline strain 179/07-11 did not hemagglutinate 
those cells despite the high HA titers obtained with mouse red 
blood cells (Table 2). Using mouse erythrocytes, RDE activity 
(loss of the HA pattern formed at 4 °C) was observed for the 
bovine and bubaline isolates (Table 2). 


Genetic characterization of the bubaline CoV 


A 9.6-kb region encompassing the entire 3’ end of the viral 
genome (from the 32-kDa to nucleocapsid protein genes) was 
determined through PCR amplification and subsequent sequenc- 
ing of overlapping fragments. At the 3’ end of viral RNA, 
bubaline CoV strain 179/07-11 had the same genomic organiza- 
tion of other ruminant CoVs (Fig. 2). Nine ORFs were identified 
by both ORF Finder program and sequence comparison with 
reference ruminant CoV sequences. All predicted ORFs but the 
4.8-kDa gene were preceded by a repeated intergenic sequence, 
CUAAAC or CCAAAC, which is predicted to interact with the 


4.9x0a 4.8k0a 


BuCoV 278 424 1363 
BCoV 278 424 1363 
GiCoV 278 424 1358 
ACoV 278 424 1363 
SACoV 278 424 1363 


4345 109 84 230 207 448 
4345 109 & 230 207 448 
2938 109 & 230 207 448 
2945 109 8& 230 207 448 
2938 109 8 230 207 448 


Fig. 2. Schematic comparison of the genomes (3’ end) of different ruminant CoVs. RdRp, RNA dependent-RNA polymerase. Below the diagram, the length in amino 
acids is reported for the encoded proteins of bubaline coronavirus 179/07-11 (BuCoV), bovine coronavirus Mebus (BCoV), giraffe coronavirus US/OH3/2003 
(GiCoV), alpaca coronavirus (ACoV) and sable antelope coronavirus US/OH1/2003 (SACoV). 
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Table 3 
Coding potential and putative transcription regulatory sequences of the 3’ end of 
the BuCoV genome 


Putative gene Putative TRS 


Gene segment Coding sequence Start nt position TRS sequence 


32 kDa 200-1036 187 CUAAAC 
HE 1048-2322 1033 CUAAAC 
S 2337-6428 2330 CUAAAC 
4.9 kDa 6418-6549 6095 CCAAAC 
4.8 kDa 6585-6722 - Not detected 
12.7 kDa 6806-7135 6725 CCAAAC 
E 7122-7376 6993 CCAAAC 
M 7390-8082 7381 CCAAAC 
N 8092-9438 8079 CUAAAC 
I 8153-8776 8079 CUAAAC 


viral polymerase during the discontinuous transcription of the 
negative strand subgenomic RNA of the Nidovirales members 
(Table 3). 

Table 4 shows the unique aa changes encountered in the 
encoded proteins of strain 179/07-11 with respect to other ru- 
minant CoVs. The predicted 32-kDa nonstructural protein (nsp) 
had the same length (278 aa) as most other bovine-like CoVs, 
including ruminant CoVs giraffe coronavirus (GiCoV), alpaca 
coronavirus (ACoV) and sable antelope coronavirus (SACoV), as 


N. Decaro et al. / Virology 370 (2008) 213-222 


well as CRCoV and HCoV-OC43. Three aa changes were found 
to be unique to strain 179/07-11. The highest aa identity (98.9%) 
of the 32-kDa protein of strain 179/07-11 was found to enteric and 
respiratory BCoVs, as well as to ACoV. The hemagglutinin- 
esterase (HE) protein (424 aa) showed three unique aa mutations 
in comparison with other ruminant CoVs and it was more closely 
related (99% aa identity) to BCoVs DB2 and ENT and SACoV. 
By analysis with the NetNglyc server, the HE of strain 179/07-11 
was found to contain nine potential N-glycosylation sites, ana- 
logously to BCoV reference strains. Moreover, the predicted site 
for neuraminidate-O-acetylesterase activity, FGDS, was detected 
at the N-terminus. The spike (S) protein was long 1363 aa and 
matched the best aa identity (98.1%) to BCoV-E-AH65 and 
SACoV-US/OH1/03. Twenty-one potential N-glycosylation 
sites were identified throughout the protein, one more than other 
ruminant CoVs, due to the presence of an additional N- 
glycosylation site, NLS, at position 1260. No deletions or insertions 
were observed in strain 179/07-11 compared with bovine-like 
CoVs. Also, the aa stretch KRRSRR, responsible for the proteolytic 
cleavage of the S protein at residue 768 into subunits S1 and S2, 
was conserved. In comparison to ruminant CoV reference strains, 
16 unique residues were identified in the S protein of strain 179/07- 
11. Nsp 4.9 kDa was 43-aa long as for BCoVs Mebus, Quebec and 
DB2 and HECV-4408, whereas in other ruminant CoVs it is 


Table 4 
Unique amino acid residues observed in structural and nonstructural proteins encoded by the 3’ end of the bubaline CoV genome in comparison with other ruminant 
CoVs 
CoV strain 32 kDa* HE Ss 

205 217 275 108 158 210 143 147 151 157 257 318 328 525 608 785 805 828 888 909 927 1260 
BuCoV-179/07-11 K N FF H P N Y F I VN T F Y G kK D TT kK R A N 
BCoV-Mebus M kK V Y S K H L L I T I L H D N _E A N kK § D 
BCoV-ENT M kK V Y S K H L L I Ab I Lb H D N _E A N KS D 
BCoV-Quebec M kK V Y S K H L L I T I L H D N _E A N kK § D 
BCoV-DB2 M K V Y S K H L L I T I L H D N _E A N KS H 
GiCoV-US/OH3/2003 M K VY S K H L L I T I L H D N _E A N KS D 
ACoV M kK V Y S K H L I I T I L H D N _E A N K § D 
SACoV-US/OH1/2003 M K VY S K H L 16, I T I 15 H D N _E A N KS D 
CRCoV-4182 M K VY S kK H L is I T I i Y G N EE VN KS D 
HECV-4408 NA NA NA Y S K H L L I Ny) I L H G N _ E A N R A OD 
HCoV-OC43 M K V Y A _ R R S Vv tT I Le D D N _ E A N KS D 
PHEV-VW572 M NP NP Y_ P kK H I L I S VeoL T D N E T N R S D 
CoV strain 4.9 kDa? 4.8 kDa? 12.7kDa M N I 

24 32 40 2 9 16 28 29 40 = 41 70 58 158 152 354 56 116 122 = 132 

BuCoV-179/07-11 P A T Q D K G F N W F L S H S K I S M 
BCoV-Mebus L E M P G M Cc S 40 Cc S I G Q N R T L R 
BCoV-ENT L NP NP P G M Vv S T Cc S I G Q N R T L R 
BCoV-Quebec L E M NP G M Cc S T Cc S I G Q N R T L R 
BCoV-DB2 L E M P G M Vv Ss so Cc S I G Q N R A L R 
GiCoV-US/OH3/2003 L NP NP P Vv M Vv S NP NP S§S I G Q N R T L R 
ACoV L NP NP P Vv M Vv S T S S I G Q N R T L. R 
SACoV-US/OH1/2003 L NP NP P Vv M Vv S NP NP §S I G Q N R T L R 
CRCoV-4182 NP NP NP NP NP NP NP NP NP NP S I G Q N R T L R 
HECV-4408 L E T Q D M G S N W S I G Q N R T LL R 
HCoV-OC43 NP NP NP NP NP NP NP NP NP NP S I G Q N NP T L R 
PHEV-VW572 F NP NP P NP NP NP NP NP NP S I G Q N R T L R 


NA, not available; NP, not present. 


* Nsp 32 kDa of HECV-4408 is not available in the GenBank database, whereas it is truncated in PHEV-VW572. 
> Nsp 4.9 kDa and 4.8 kDa are not present in the CRCoV-4182 and HCoV-OC43 genomes. Nsp 4.9 kDa is truncated in PHEV-VW572. 
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truncated in the C-terminus (only 29 aa in length). The length of nsp 
4.8 kDa was 45 aa as for HECV-4408 and ruminant CoVs with the 
exception of GiCoV and SACoV (38 aa). Nsp 12.7 kDa showed the 
same size as for HECV-4408 and ruminant CoVs (109 aa), 
although BCoV Quebec exhibits a truncate form of the protein (87 
aa). In the three nonstructural proteins, strain 179/07-11 displayed 
the highest identity to HECV-4408, with three, seven and one 
unique residues in nsp 4.9 kDa, 4.8 kDa and 12.7 kDa, respectively. 
The envelope (E) protein was 84-aa long exactly as BCoV and 
GiCoV but 1-aa shorter and 2-aa longer than the analogous protein 
of ACoV and SACOoV, respectively. No residue unique to the 
bubaline strain was detected at this level. A 100% aa identity in the 


A 


E protein was found between the bubaline strain and GiCoV and 
most BCoV strains, with the exception of BCoV-Mebus (98.8%), 
whereas the bubaline isolate was slightly less related to ACoV and 
SACoV. The membrane (M) protein was 220-aa long and contained 
two unique residues, with the I58L change shared with HECV- 
4408. In this protein the highest identity (99.5%) was matched to 
HECV-4408. One potential N-glycosylation and four potential O- 
glycosylation sites were detected at the N-terminus of the M 
protein, as for all bovine-like CoVs including the closely related 
HECV-4408, with which the bubaline isolates shared an aa 
mutation with respect to other ruminant CoVs. However, both 
programs indicated that M gene sequence of the ruminant CoVs 


137 ACoV 
BCoV-ENT; GiCoV; SACoV 
BCoV-R-AH65; BCoV-R-AH65-TC 


32! BCoV-DB2 


BCoV-LUN 


35] BCoV-E-AH65 
4a BCOV-E-AH187 


BCoV-E-AH65-TC 


20} HECoV-4408 


BCoV-R-AH187 
BCoV-Quebec 


54 | 6! 'BCoV-Mebus 


100 


———1 
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BuCoV-179/07-11 
CRCoV-4182 
HCoV-0C43 


PHEV-VW572 
MHV-A59 
SDAV-681 
HCoV-HKU1 


BuCoV-179/07-11 
HECoV-4408 


58 |' BCoV-DB2 


BCoV-E-AH187 
BCoV-ENT 


441 BCoV-LUN; E-AH65-TC; R-AH65-TC; R-AH187; ACoV 
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PHEV-VW572 


HCoV-HKU1 


-——_ 
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Fig. 3. Neighbor-joining trees based on the hemagglutinin-esterase (A) and membrane (B) proteins of group 2 CoVs. For phylogenetic tree construction, the following 
CoV strains were used (GenBank accession numbers are reported in parentheses): BCoV-Mebus (U00735), Quebec (AF220295), DB2 (DQ811784), ENT 
(AF391541), LUN (AF391542), E-AH65 (EF424615), R-AH65 (EF424617), E-AH65-TC (EF424616), R-AH65-TC (EF424618), E-AH187 (EF424619), R-AH187 
(EF424620); GiCoV-US/OH3/2003 (EF424623); ACoV (DQ915164); SACoV-US/OH 1/2003 (EF424621); CRCoV-4182 (DQ682406); HCoV-OC43 (NC_005147); 
HECV-4408 (HE, L07747; M, AY316299); PHEV-VW572 (DQ011855); MHV-A59 (AY700211); SDAV (AF207551); HCoV-HKU1 (NC_006577). A statistical 
support was provided by bootstrapping over 1000 replicates. The scale bars indicate the estimated numbers of amino acid substitutions per site. 
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might not contain a signal peptide. Proteins without signal pep- 
tides are unlikely to be exposed to the O-glycosylation and N- 
glycosylation machineries and thus may not be glycosylated in vivo 
even though they contain potential motifs. The nucleocapsid (N) 
protein of strain 179/07-11 had a length of 448 aa and was closely 
related (99.3% of aa identity) to BCoV-DB2. Two unique aa 
changes were found in comparison to other bovine-like viruses. The 
I protein (207 aa) encoded by the N-gene internal ORF showed four 
unique residues and the best aa identity (96.1%) to BCoV reference 
strains Mebus, Quebec and DB2. When the 9.6-kb sequence 
corresponding to the 3’ end of CoV genome was analyzed, the 
bubaline strain showed over 97.6% nt identity to BCoV strains and 
other ruminant CoVs, whereas it was less genetically related to 
other bovine-like CoVs, including CRCoV (96.5% nt identity), 
HCoV-0C43 (91.4%), HECV-4408 (85.1%) and PHEV-VW572 
(83.1%). 

Phylogenetic analysis with the structural protein HE showed 
that the bubaline strain 179/07-11 clusters together with the 
bovine-like CoVs, being more related to ruminant viruses 
(Fig. 3). Similar phylogenetic trees were generated for the other 
proteins (data not shown) but the M protein, where strain 179/ 
07-11 forms a separate cluster with HECV-4408 into the bovine 
subgroup (Fig. 3). 


Evaluation of the genetic distance among ruminant Covs 


In order to better define the taxonomic position of the bubaline 
CoV 179/07-11, a sequence comparison between BCoV reference 
strain ENT and ruminant CoVs GiCoV, ACoV, SACoV and the 
bubaline strain was conduced in the nonstructural and structural 
proteins (Table 5). For such a comparison, strain ENT was chosen 
among the available BCoV sequences as it is an enteric strain 
isolated very recently (Chouljenko et al., 2001). The results 
showed that strain 179/07-11 is less related to BCoV-ENT than to 
other ruminant CoVs in all major structural and nonstructural 
proteins, whereas closer relationship (98.8% aa identity) to 
BCoV-ENT with respect to SACoV (97.6% aa identity) was 
found only in the E protein. Analysis of the nt sequence cor- 


Table 5 

Sequence identity (%) of ruminant CoVs to reference BCoV-ENT in 
nonstructural and structural proteins and in the full-length sequence of the 
genomic RNA 3’ end (9.6 kb) 


Gene CoV strain 
segment BuCoV GiCoV ACoV  SACoV 
179/07-11  US/OH3/2003 US/OH1/2003 
32 kDa® 98.2 98.5 99.2 98.5 
HE? 99.0 100 99.7 100 
s? 98.0 98.8 99.1 99.3 
4.9 kDa® 55.8 100 100 100 
4.8 kDa® 75.5 82.2 93.3 82.2 
12.7 kDa® 98.1 99.0 99.0 98.1 
EA 98.8 100 98.8 97.6 
M’ 98.2 99.5 99.5 99.5 
N’ 98.4 99.5 99.7 99.5 
I 95.1 98.0 99.0 98.5 
9.6-kb 3’ end’ «97.8 99.1 99.3 99.2 


* Amino acid identity. 
> Nucleotide identity. 


responding to the entire 3’ end of CoV genome confirmed that the 
bubaline strain has a lower genetic relationship (97.8% nt identity) 
to the reference BCoV than other recently isolated ruminant CoVs 
(99.1—99.3% nt identity). Confirmation of the conclusions drawn 
from Table 5 was derived with the use of SimPlot analysis. 


Discussion 


In the present report, a bovine-like CoV in water buffalo (B. 
bubalis) has been described, along with its biological and 
genetic characterization. The bubaline strain was detected in the 
feces of buffalo calves with severe gastroenteritis, causing the 
death of at least 30 animals. A virus (strain 179/07-11) anti- 
genically related to BCoV was isolated on HRT-18 cells, 
producing CPE typical of BCoV and being recognized by a 
BCoV-specific serum in an IF assay. Analysis of the biological 
properties of the bubaline strain showed that, unlike classical 
BCoV strains (Fukutomi et al., 1999), the virus was not able to 
replicate efficiently on MDBK cells or cause hemagglutination 
with chicken erythrocytes. However, Benfield and Saif (1990) 
demonstrated that winter dysentery strains of BCoV from adult 
cattle were also unable to replicate in MDBK cells. Among wild 
ruminant CoVs, only few strains have been characterized at a 
biological level. GiCoV displayed HA activity also with chick- 
en erythrocytes at a temperature of 4 °C, whereas no data 
are available about its replication efficiency on MDBK cells 
(Hasoksuz et al., 2007). Using chicken erythrocytes, other 
ruminant CoVs, including the waterbuck, sambar deer and 
white-tailed deer isolates and some BCoV strains, exhibited 
very low HA titers at 4 °C and no HA activity at 37 °C. Thus, 
the loss of ability to replicate on MDBK cells and to agglutinate 
chicken erythrocytes appears as characteristic of the bubaline 
CoV with respect to the other ruminant CoVs. 

Sequence analysis of the 3’ end of the viral genome showed 
that strain 179/07-11 has a genomic organization similar to 
BCoV, including the presence of three ORFs encoding for small 
nonstructural proteins between the S and E genes. Moreover, the 
virus has a close genetic relatedness to the bovine CoV sub- 
group in the major structural and nonstructural proteins. How- 
ever, as shown by sequence and SimPlot analyses, the genetic 
distance between the bubaline virus and reference BCoV strain 
ENT was approximately the same or even higher than observed 
between the BCoV strain and other ruminant CoVs. Unique aa 
substitutions were detected in all proteins of the bubaline virus 
but they were mainly within the S protein (Table 4). As ob- 
served in BCoV (Chouljenko et al., 1998), GiCoV (Hasoksuz 
et al., 2007) and ACoV (Jin et al., 2007), these changes occurred 
mostly in the S1 subunit which has been proven to be involved 
in receptor binding and host specificity in the well-characterized 
MHV (de Haan et al., 2006). Certain residues in the S protein 
had been associated to enteric or respiratory tropism of BCoV, 
including residues 531, 769 and 1026 (Chouljenko et al., 1998). 
The bubaline strain displayed the same aa residues encountered 
in enteric BCoVs in only one position (D531). The remaining 
two aa residues (S769, G1026) suggest a potential respiratory 
tropism, which is in contrast with the exclusive enteric clini- 
cal pattern observed in water buffalo. However, more recent 
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analysis of two enteric and respiratory BCoV pairs failed to 
identify specific aa residue determinants in the S proteins as 
well as in other structural and nonstructural proteins, except for 
a T22I mutation in nsp 4.8 kDa (Zhang et al., 2007). In- 
terestingly, strain 179/07-11 showed the same aa determinant 
T22 as the enteric BCoV strains. Other authors have identified a 
truncated form of nsp 4.9 kDa (29 instead of 43 aa) as a genetic 
marker for the respiratory tropism (Gelinas et al., 2001; Vijgen 
et al., 2006). The bubaline strain presented an intact nsp 4.9 kDa 
as enteric BCoVs, but even this aa deletion has been recently 
shown not to be involved in tissue tropism (Zhang et al., 1994). 

Based on the unique biological properties and the more distant 
relatedness of strain 179/07-11 to BCoV with respect to other 
ruminant CoVs, we propose to designate this strain as prototype 
of a host-range variant of BCoV, namely bubaline CoV (BuCoV). 
With regard to the its possible origin, BuCoV may have arisen 
through interspecies transmission of a BCoV strain from cattle to 
water buffaloes, although no recent contacts with cattle or other 
ruminants were reported. This hypothesis is supported by the high 
genetic relatedness between BuCoV and BCoV. Moreover, it 
should be considered that a bovine origin has been strongly 
suggested for other group 2 CoVs less genetically related to 
BCoV, such as HCoV-OC43 (Vijgen et al., 2005), PHEV (Vijgen 
et al., 2006) and CRCoV (Erles et al., 2007). 

Accumulation of point mutations, as well as small insertions 
and deletions in coding and non-coding sequences, are the 
dominant forces in the microevolution of plus-stranded (+) RNA 
viruses, resulting in proliferation of virus strains, serotypes and 
subtypes (Dolja and Carrington, 1992). Coronaviruses are thought 
to mutate at high frequency like most RNA viruses as a con- 
sequence of high error rates of the RNA polymerase that are 
predicted to accumulate several base substitutions per round of 
replication (Jarvis and Kirkegaard, 1991; Lai and Holmes, 2001). 
Changes in tissue tropisms and/or interspecies transmission of 
CoVs occur through genetic variations in structural and/or non- 
structural proteins (Laude et al., 1993; Vennema et al., 1998; Guan 
et al., 2003; Rottier et al., 2005; Song et al., 2005; Vijgen et al., 
2005; Decaro et al., 2007c). Genetic determinants that may have 
caused the cattle-to-buffalo interspecies transmission have not 
been identified so far. Although this study has detected a CoV 
strain in buffalo calves with severe diarrhea, the pathogenicity of 
this virus and its etiologic role in enteric disease of water buffalo 
still have to be determined more definitively. In addition, epide- 
miological studies could assess whether BuCoV is widespread 
among water buffalo herds as well as whether cross-species 
transmission between buffalo and cattle occurs mainly in areas 
where both closely related ruminant species are raised intensively. 


Materials and methods 
(RT-)PCR assay for screening for bovine viral pathogens 


Nucleic acids for (RT-)PCR assays were purified using the 
commercial kits DNeasy Tissue Kit (QIAGEN S.p.A., Milan, 
Italy) and QIAamp® Viral RNA Mini Kit (Qiagen S.p.A., Milan, 
Italy) from the fecal samples and intestinal contents and QIAamp® 
RNeasy Mini Kit (Qiagen S.p.A.) from the tissue samples. 


RNA extracts were used for detection of BCoV (Erles et al., 
2003), toroviruses (Hoet et al., 2002), rotaviruses (Gouvea et al., 
1994), caliciviruses (Jiang et al., 1999), bovine viral diarrhea 
virus (BVDV) (Sullivan and Akkina, 1995) and bovine 
respiratory syncytial virus (Valarcher et al., 1999). Detection 
of bovine herpesvirus types 1 (BoHV-1) (Vilcek, 1993) and 4 
(BoHV-4) (Boerner et al., 1999) was carried out on the DNA 
templates. RT-PCR and PCR assays were performed using 
SuperScript™ One-Step RT-PCR for Long Templates (Life 
Technologies) and LA PCR Kit Ver. 2.1 (TaKaRa Bio Inc., 
Shiga, Japan), respectively. 


Bacteriological and parasitological investigations 


The samples were examined for bacterial and parasitic 
pathogens by standardized methods. For bacteriological in- 
vestigations, the fecal samples were plated out on 5% sheep 
blood agar and cultured aerobically at 37 °C for 24 h to exclude the 
presence of aerobic pathogens. One hundred milligrams of fecal 
samples was resuspended in 900 l of fluid thioglycolate me- 
dium (FTG). Ten-fold dilutions (10 * to 10°) were subsequently 
plated onto 5% sheep blood agar and egg yolk agar with D- 
cycloserine 400 t1g/ml. Bacteria were allowed to grow overnight 
at 37 °C in anaerobic condition. Detection of the most common 
enteric parasites was achieved using zinc sulfate flotation. The 
Ziehl Nielsen staining was also performed on fecal samples or 
intestinal sections for detection of Cryptosporidium spp. 


Real-time RT-PCR for quantification of bovine-like CoVs 


A real-time RT-PCR assay based on TaqMan technology 
(G. Elia et al., manuscript in preparation) was used to quantify 
the viral load in samples tested positive for bovine-like CoVs by 
conventional RT-PCR (Erles et al., 2003). The TaqMan assay 
had been proven to be able to detect BCoV as well as the closely 
related CRCoV (Elia et al., manuscript in preparation). BCoV-like 
RNA copy numbers were calculated on the basis of the standard 
curves generated by 10-fold dilutions ofa synthetic RNA obtained 
by in vitro transcription of a plasmid containing the M gene of 
BCoV strain 339/07 (Decaro et al., 2007b). Reverse transcription 
was carried out using GeneAmp® RNA PCR (Applied Biosys- 
tems, Applera Italia, Monza, Italy), following the manufacturer’s 
recommendations. The quantitative assay targeting the M gene 
was conducted in a 50-ul reaction mixture containing 25 ul of 
IQ™ Supermix (Bio-Rad Laboratories Srl), 600 nM of primers 
BCoV-F (CCTTCATATCTATACACATCAAGTTGTT) and 
BCoV-R (ACCAGCCATTTTAAATCCTTCA), 200 nM of 
probe BCoV-Pb (6FAM CCTTCATATCTATACACAT- 
CAAGTTGTT-BHQ]1) and 20 tl of c-DNA. The thermal profile 
consisted of activation of iTfaq DNA polymerase at 95 °C for 
10 min followed by 45 cycles of denaturation at 95 °C for 15 s and 
annealing extension at 60 °C for | min. 


Virus isolation 


Virus isolation was carried out on the buffalo fecal sample 
containing the highest RNA titer of BCoV-like coronavirus 
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(strain 179/07-11), as determined by real-time RT-PCR. Feces 
were homogenized in Dulbecco’s minimal essential medium 
(D-MEM) containing antibiotics (penicillin 5000 IU/ml, strep- 
tomycin 2500 pg/ml, amphotericin B 10 pg/ml). After cen- 
trifugation at 3000xg for 15 min, the supernatant was used to 
inoculate confluent monolayers of HRT-18 and MDBK cells in 
the presence of trypsin (5 ug/ml). Viral growth was monitored 
by an immunofluorescence (IF) assay using a BCoV-positive 
bovine serum and a rabbit anti-bovine IgG conjugated with 
fluorescein isothiocyanate (Sigma Aldrich srl, Milan, Italy). 


HA and RDE activity 


Two-fold dilutions (starting from dilution 1:2) of the super- 
natant of the CoV-infected HRT-18 cells tested positive by the IF 
assay were made in phosphate buffered saline (PBS, pH 7.2) using 
96-well V-plates. Mouse (0.8%) or chicken (0.4%) erythrocytes 
were added to each dilution. Results were read after 1 h at 4 °C or 
37 °C and expressed as the reciprocal of the highest dilution 
producing HA. The plates were then incubated for additional 4 h at 
37 °C to assess RDE activity, which was expressed as the 
reciprocal of the highest dilution of virus resulting in complete 
disappearance of HA (Hasoksuz et al., 1999). 


PCR amplifications of the 3' end of the BCoV-like virus 


The 3’ end of the genome of strain 179/07-11 was amplified 
from the original fecal sample using SuperScript™ One-Step 
RT-PCR for Long Templates (Life Technologies, Invitrogen. 


Table 6 


Primers used for RT-PCR amplification and sequence analysis 


Milan, Italy), according to the manufacturer’s instructions. 
Thirteen partially overlapping fragments encompassing from the 
very 3’ end of ORF1b to the 5’ end of the untranslated region 
(UTR) were amplified using primer pairs designed on conserved 
regions among BCoV-like group 2 coronaviruses (Table 6). 


Sequence analysis and phylogeny 


The PCR-amplified products were sequenced by Genome 
Express (Meylan, France) and the obtained sequences were 
assembled and analyzed using the BioEdit software package 
(Hall, 1999) and the NCBI’s (htttp://www.ncbi.nlm.nih.gov) 
and EMBL’s (http://www.ebi.ac.uk) analysis tools. Additional 
RT-PCR assays and sequencing attempts were performed to 
close gaps between assembled contigs using strain-specific 
primers. ORFs contained in the amplified genomic region were 
determined either with the ORF Finder tool of NCBI or on the 
basis of the similarity to known coronavirus proteins. The ORFs 
identified in this manner were translated and the predicted 
amino acid (aa) sequences were saved as individual files for 
further analyses. Putative N-glycosylation sites were predicted 
using the NetNglyc 1.0 Server program (http://www.cbs.dtu.dk/ 
services/NetNGlyc/), whereas the NetOGlyc 3.1 Server program 
(http://www.cbs.dtu.dk/services/NetOGlyc/) was used for predic- 
tion of O-linked glycosylation sites. Phylogenetic and molecular 
evolutionary analyses were conducted using Mega3 (Kumar et al., 
2004). Phylogenetic trees, based on the nonstructural 32-kDa 
protein, on the structural spike (S), envelope (E), membrane (M) 
and nucleocapsid (N) proteins and on the genomic RNA 3’ end 


Primer Sequence 5’ to 3’ Sense Position * Amplicon size (bp) 
BCV-22001F TAGACTTGAAATAGTTAAGCTTGGTG + 22,001—22,026 893 
BCV-22894R AAATTAGCTTCACGAGCTATATATGC = 22,868—22,893 

BCV-22769F TATCGCAGCCTTACTTTTGTTAATG a 22,769—22,793 527 
BCV-23295R CGAAAATAACAGTACGGGGGTTGACA - 23,270—23,295 

BCV-23112F TACCCTCTGGTAATTATTTAGCCATTTCA ate 23,112—23,140 776 
BCV-23887R TTCCCTTCAGTGCCATATTACGATATGT = 23,860—23,887 

BCV-23510F TATGATCCGCTACCAATTATTTTGCTTGGCA + 23,510—23,540 817 
BCV-24326R ACAACACCAGTGTCTGTAAAATATGCA = 24,300—24,326 

BCV-24182F TAGAACTATGGCATTGGGATACAGGTGTTG t 24,182—24,211 1254 
BCV-25435R TACACCTATCCCCTTGTAAACAAGAGTA = 25,409—25,435 

BCV-25301F ACTTAGTTGGCATAGGTGAGCACTGTTC + 25,301—25,328 851 
BCV-26151R ACATGCTACATAATCACCACAGACAA - 26,126—26,151 

BCV-26028F TTTGTATGAAATTCAAATACCTTCAGAG te 26,028—26,055 877 
BCV-26904R GTCTATCTGAGCTTGCGCTTCAAGAGCA = 26,877—26,904 

BCV-26760F TAAAATTCAAGCTGTTGTTAATGCAAAT + 26,760—26,787 1059 
BCV-27818R GCTCGACCTAAATGGGTCTTATAATTAGA = 27,791—27,818 

BCV-27667F GGTGGTTGTTGTGATGATTATACTGGACA Sa 27,667—27,695 853 
BCV-28519R ACTACAACTATTATAACCAATAAACACAT = 28,49 1—28,519 

BCV-28380F TCTTAGCTGTTGACTTTATTACCTGG + 28,380—28,405 919 
BCV-29298R ACATAAACAGCAAAACCACTAGTATCGCC - 29,270—29,298 

BCV-29160F AAGGTATAAAACTAGGTACTGGCTAT ote 29,160—29,185 876 
BCV-30035R TGCGCGATCCTGCACTAGAGGCTCTAC = 30,009—30,035 

BCV-29900F GAGGCTATTCCGACTAGGTTTCCGCCT + 29,900—29,926 827 
BCV-30726R GTGTCTTCAGTAAAGGGCTCATCCATC = 30,700—30,726 

BCV-30610F TGATAATATAAGTGTTGCAGCGCCCAAA =f 30,610—30,637 419 
BCV-31028R GTGATTCTTCCAATTGGCCATAATT = 31,004—31,028 


* Primers position is referred to the sequence of BCoV strain ENT (accession number AF391541). 
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(9.6 kb) of CoV strain 179/07-11, were elaborated using both 
parsimony and neighbor-joining methods, supplying a statistical 
support with bootstrapping over 1000 replicates. 


SimPlot analysis 


The SimPlot program (version 3.2) was used to analyze the 
genetic distance of the 3’ end of the genomes of ruminant CoVs, 
including the bubaline strain, in reference to the same region of 
BCoV strain ENT and this genetic distance was plotted versus 
nucleotide (nt) positions (Lole et al., 1999). 


Nucleotide sequence accession number 


The nt sequence of the bubaline strain 179/07-11 has been 
deposited in GenBank under accession number EU019216. 
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